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ABSTRACT 

A  high-resolution,  multi-level,  primitive  equation  ocean  model  is  used  to 
investigate  the  roles  of  wind  forcing  and  irregular  coastline  geometry  in  the  generation  of 
currents,  eddies,  jets  and  filaments  in  the  Canary  Current  System  (CCS)  from  30°  N  to 
42.5°  N.  To  study  the  generation,  evolution,  and  sustainment  of  the  currents,  eddies,  jets 
and  filaments  in  the  CCS,  the  model  is  forced  from  rest  using  seasonal  climatological 
winds  and  a  realistic  coastline.  Results  of  the  experiment  show  that  wind  forcing  alone  is 
capable  of  generating  surface  currents,  undercurrents,  meanders,  eddies,  and  filaments. 
Preferred  eddy  generation  locations,  enhanced  growth  of  meanders,  eddies,  and  filaments 
are  seen.  The  features  produced  by  the  model  are  consistent  with  available  observations 
of  the  CCS. 


VI 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 1 

II  MODEL  DESCRIPTION 5 

A.  MODEL  EQUATIONS 5 

B.  TYPE  OF  WIND  FORCING 7 

C.  EXPERIMENTAL  DESIGN 8 

IE.        RESULTS  FROM  THE  MODEL  SIMULATION 11 

A.  SPIN-UP  PHASE 11 

B.  QUASI-EQUILBRIUM  PHASE 12 

C.  COMPARISONS  OF  MODEL  RESULTS  WITH  OBSERVATIONS..  15 

1.  Comparison  of  Ocean  Currents 16 

2.  Comparison  of  Eddies 17 

3.  Comparison  of  Upwelling 17 

4.  Comparison  ofFilaments 18 

IV.        SUMMARY 21 

APPENDIX.  METHOD  OF  SOLUTION 67 

LIST  OF  REFERENCES 73 

INITIAL  DISTRIBUTION  LIST 79 


vn 


Vlll 


LIST  OF  FIGURES 

1.  The  model  domain  for  the  Canary  Current  System  (CCS)  is  bounded  by  30°  N  to 
42.5°  N,  5°  W  to  17.5°  W.  The  cross-shore  (alongshore)  resolution  is  9  km  (1 1  km). 
Geographic  locations  and  prominent  features  are  labeled 25 

2.  Climatological  (1980-1989)  ECMWF  winds  in  m/s  for:  (a)  February,  (b)  July,  (c) 
September,  and  (d)  November.  Maximum  wind  vector  is  10  m/s 26 

3.  Figure  3.  Temperature  contours  and  velocity  vectors  at  30  m  depth  at  days  (a)  150, 
(b)  195,  (c)  285  and  (d)  360.  The  contour  interval  is  1°C.  To  avoid  clutter,  the 
velocity  vectors  are  plotted  every  third  grid  point  in  the  cross-shore  direction  and 
every  fourth  grid  point  in  the  along-shore  direction.  Maximum  current  velocity  is  100 
cm/s 30 

4.  Cross-shore  sections  of  isotherms  on  day  195  at  (a)  33°  N  and  (b)  39°  N.  Contour 
interval  is  1°  C 34 

5.  Cross-shore  sections  of  meridional  velocity  (v)  on  day  240  at  (a)  39.5  °  N,  and  (b) 
35.5  °  N.  Solid  lines  indicate  poleward  flow,  while  dashed  lines  indicate  equatorward 
flow.  The  contour  interval  is  2.5  cm/s  for  poleward  flow  and  5  cm/s  for  equatorward 
flow 36 

6.  Figure  6.  Temperature  contours  (b-e)  and  velocity  vectors  (a-e)  at  30  m  depth  in  the 
third  year  of  the  model  simulation  time-averaged  over  the  months  of  (a)  January,  (b) 
June,  (c)  August,  (d)  September  and  (e)  November.  The  contour  interval  is  1°  C. 
Maximum  current  velocity  is  100  cm/s 38 

7.  Cross- shore  sections  of  meridional  velocity  (v)  in  the  third  year  of  model  simulation 
time-averaged  for  (a)  January  at  39.4°  N,  (b)  June  at  38°  N  and  (c)  38.6°  N,  (d) 
August  at  37.5°  N,  and  (e)  December  at  41.5°  N.  Solid  lines  indicate  poleward  flow, 
while  dashed  lines  indicate  equatorward  flow.  The  contour  interval  is  2.5  cm/s  for 
poleward  flow  in  (a),  (d)  and  (e),  5  cm/s  for  poleward  flow  in  (b)  and  (c),  and  5  cm/s 
for  equatorward  flow 43 


IX 


8.  Time-averaged  plots  for  the  peak  upwelling  season  (July- September)  of  temperature 
contours  and  velocity  vectors  (a)  at  30  m  depth  and  (d)  at  1226  m  depth,  and  a  cross- 
shore  section  of  velocity  (v)  at  (b)  37.5°  N,  (c)  39°  N  and  (e)  37.4°  N.  The  contour 
interval  is  1°  C  in  (a).  The  maximum  current  velocity  is  100  cm/s  in  (a)  and  40  cm/s 
in  (d).  In  (b),  (c)  and  (d)  solid  lines  indicate  poleward  flow,  while  dashed  lines 
indicate  equatorward  flow.  The  contour  interval  for  poleward  flow  is  2.5  cm/s  in  (b) 
and  (c)  and  5  cm/s  in  (e).  The  contour  interval  for  equatorward  flow  is  5  cm/s 48 

9.  (a)  Mean  and  (b)  eddy  kinetic  energy  at  30  m  depth  for  the  peak  upwelling  season 
(July  through  September).  Contour  interval  is  100  cm2/s2 53 

10.  Computer  representation  of  AVHRR  NOAA7  Sea  Surface  Temperatures  (SST)  with 
observed  eddy  formations  from  5  August  1982.  "After  Fiuza  1984." 55 

11.  Temperature  contours  and  velocity  vectors  at  30  m  showing  upwelling  in  the  third 
year  at  day  (a)  137,  (b)  152,  (c)  173,  (d)  209,  (e)  263,  (f)  290,  and  (g)  317.  The 
contour  interval  is  1  °  C.  The  maximum  current  velocity  is  1 00  cm/s 56 

12.  AVHRR  sea  surface  temperatures  (SST)  from  (a)  21  June  1995  showing  increased 
cold  SST  southwest  of  Cabo  de  Sao  Vincent  caused  by  intensified  upwelling  and 
advection  of  upwelled  water  around  the  cape  and  (b)  27  July  1 993  showing  strong 
upwelling  and  numerous  filaments  off  the  west  coast  of  the  IP.  In  (c)  SST  show 
intense  upwelling  north  of  Cabo  da  Roca  on  29  July  1977.  Figures  (a)  and  (b)  are 
from  the  Remote  Sensing  Data  Analysis  Service  (RSDAS),  Plymouth  Marine 
Laboratory,  UK.  Figure  (c)  is  from  Fiuza  (1982) 63 


LIST  OF  TABLES 

1.     Values  of  Constants  Used  in  the  Model 71 


XI 


Xll 


ACKNOWLEDGEMENT 

The  professional  guidance  and  considerable  knowledge  of  Dr.  Mary  Batteen  made 
this  study  possible.  Hugh  thanks  to  Pete  Braccio,  Mike  Cook  and  Phaedra  Green  for 
FORTRAN,  FERRET  and  MATLAB  assistance.  Thanks  also  to  NERC  Remote  Sensing 
Data  Analysis  Service,  Plymouth  Marine  Laboratory,  UK  for  the  use  of  their  high  quality 
AVHRR  SST  images  (http://www.ac.uk/rsdas). 

Finally,  I  would  like  to  thank  my  wife  for  her  patience  and  support  during  the  last 
twenty-one  months. 


xni 


XIV 


I.  INTRODUCTION 

The  Canary  Current  System  (CCS)  is  a  classical  eastern  boundary  current  (EBC) 
system  located  off  the  west  coasts  of  Northwest  Africa  and  the  Iberian  Peninsula  (IP).  It 
extends  from  -10°  N  to  45°  N,  and  forms  the  closing,  eastern  boundary  of  the  North 
Atlantic  gyre.  The  climatological  mean  CCS  consists  of  several  large-scale  currents  of 
which  the  Canary  Current  (CC)  is  predominant.  Like  other  EBCs,  the  CC  is  a  broad 
(-1000  km),  relatively  slow  (-10-30  cm/s),  equatorward,  yearlong  surface  flow  extending 
to  depths  of-  500  m  (Wooster  et  al.,  1976).  The  portion  of  the  CC  which  lies  off  the 
west  coast  of  the  IP  is  occasionally  referred  to  as  the  Portugal  Current  (e.g.,  Tomczak  and 
Godfrey,  1994).  Embedded  in  the  CC  are  narrower,  somewhat  swifter  flows  occasionally 
reaching  speeds  of  100  cm/s  (DMAH/TC,  1988)). 

Beneath  the  CC  and  near  the  coast,  a  narrower  (-10-40  km)  and  weaker  (-2-10 
cm/s)  poleward  undercurrent  is  present.  While  the  undercurrent  is  usually  strongest 
between  100  and  600  m  depth  (Haynes  and  Barton,  1990),  the  depth  and  strength  of  the 
undercurrent  can  vary  both  seasonally  and  latitudinally.  During  the  winter  the 
undercurrent  shoals  to  the  north,  occasionally  reaching  the  surface,  to  form  a  third  flow 
component  referred  to  as  the  Iberian  Current  (IC)  (Haynes  and  Barton,  1990).  The  IC  is  a 
narrow  (25-40  km),  relatively  weak  (-20-30  cm/s)  seasonal  surface  poleward  current 
found  trapped  near  the  coast  against  the  shelf  break  (Fiuza,  1980;  Frouin  et  al.,  1990; 
Haynes  and  Barton,  1990).  Geographically,  the  IC  is  normally  found  north  of  Cabo  da 
Roca  (see  Figure  1  for  geographical  locations)  but  occasionally  extends  south  nearly  to 
Cabo  de  Sao  Vincent  during  the  winter  months.  A  second  deeper  undercurrent  also  flows 
poleward  off  the  west  coast  of  the  IP.  Attributed  to  the  Mediterranean  Outflow  (MO),  this 
flow  has  at  least  two  main  cores,  a  shallow  core  at  depths  of  600-900  m  and  a  deeper  core 
at  1100-1200  m  (Ambar,  1980).  An  additional  third,  and  shallower,  poleward  core  of 
Mediterranean  water  has  also  been  shown  to  exist  near  400  m  (Fiuza,  1979;  Ambar, 
1980). 


The  CCS  is  influenced  predominantly  by  equatorward,  upwelling  favorable  winds 
produced  by  the  Azores  High.  The  Azores  High  is  a  semi-permanent  subtropical  high 
pressure  system  similar  in  nature  and  behavior  to  its  counterpart  off  the  west  coast  of 
North  America,  the  North  Pacific  Subtropical  High,  as  described  in  Nelson  (1977).  The 
center  of  the  Azores  High  migrates  meridionally  with  the  seasons,  reaching  -27°  N,  its 
southernmost  extent,  in  March  and  ridging  north  to  -33°  N  by  August.  As  a  result  of  this 
migration,  maximum  wind  stress  values  vary  temporally  at  given  locations.  The  summer 
mean  east-west  pressure  contrast  between  Portugal  and  the  center  of  the  Azores  High  is 
~8  mb.  During  the  winter,  the  contrast  lessens  to  ~1  mb.  This  pressure  difference  results 
in  correspondingly  strong  northerly  to  northwesterly  winds  during  the  summer  and  weaker 
northwesterly  to  even  slight  southerly  winds  off  of  northern  Portugal  and  northwest  Spain 
during  the  winter.  The  shift  of  maximum  wind  stress  also  causes  the  upwelling  favorable 
winds  to  shift  from  -27°  N  near  the  Canary  Islands  in  January,  to  -43°  N  off  the  IP  by 
July.  (Fiuza,  1982) 

Recent  observations  have  shown  highly  energetic,  mesoscale  features  such  as  jet- 
like  surface  currents,  meanders,  eddies,  and  filaments  superimposed  on  the  broad, 
climatological  mean  flow  of  the  CC  and  other  EBCs.  Specifically,  satellite  images  of  sea 
surface  temperatures  have  shown  several  filaments  extending  off  the  coast  of  Portugal  and 
northwest  Spain  (Fiuza  et  al.,  1989)  as  well  as  off  Cape  Ghir  in  northwest  Africa  (Van 
Camp  et  al.,  1991;  Hagen  et  al.,  1996).  Observations  have  also  shown  pairs  of 
anticyclonic  and  cyclonic  mesoscale  eddies  on  the  order  of  100  km  off  the  coast  of  Iberia 
(Fiuza,  1984).  Observed  during  periods  of  predominantly  equatorward,  upwelling 
favorable  winds,  these  features  seem  to  be  located  near  prominent  coastal  features  such  as 
capes.  These  observations  provide  evidence  that  wind  forcing  and  coastal  irregularities 
could  be  important  mechanisms  in  the  formation  and  sustainment  of  many  of  the 
mesoscale  features  found  in  EBC  regions. 

Over  the  past  few  decades  numerous  wind  forcing  models  of  EBCs,  particularly  of 
the  California  Current  System,  have  been  conducted.    Early  work  included  steady  wind 


stress  (Pedlosky,  1974)  and  transient  wind  forcing  (Philander  and  Yoon,  1982).  The 
response  of  reduced  gravity  models  to  realistic  coastal  winds  was  investigated  by  Carton 
(1984)  and  Carton  and  Philander  (1984).  McCreary  et  al.  (1987)  conducted  a  series  of 
experiments  using  a  linear  model  with  both  transient  and  steady  wind  forcing  in  the 
California  Current  System.  In  each  of  these  previous  models,  relatively  weak  currents  (5- 
10  cm/s)  were  generated  and  no  eddies  or  filaments  developed. 

Recent  modeling  studies  have  focused  on  the  driving  mechanisms  of  complex 
features,  such  as  filaments,  observed  in  EBC  regions.  Ikeda  et  al.  (1984a,  b)  and 
Haidvogel  et  al.  (1991)  studied  baroclinic  and  barotropic  instability,  coastal  irregularities, 
and  bottom  topography  as  possible  mechanisms,  while  Batteen  et  al.  (1989),  McCreary  et 
al.  (1991),  Pares-Sierra  et  al.  (1993),  and  Batteen  (1997)  studied  wind  forcing  as  a 
possible  generative  mechanism.  Recently,  a  multilevel  primitive  equation  (PE)  model  of 
the  California  Current  System  from  35°  N  to  47.5°  N  was  used  by  Batteen  (1997)  to  study 
the  effects  of  seasonal  winds  and  coastal  irregularities  and  by  Batteen  and  Vance  (1998)  to 
study  the  additional  effects  of  thermohaline  gradients.  In  a  more  recent  study  Batteen  and 
Monroe  (1998)  studied  the  contribution  of  seasonal  wind  forcing,  thermohaline  gradients 
and  irregular  coastline  geometry  to  the  generation  of  eddies  and  filaments  for  the  entire 
California  Current  System,  i.e.,  from  22.5°  N  to  47.5°  N.  The  results  showed  that  wind 
forcing  was  the  dominant  process  responsible  for  many  of  the  observed  features  of  the 
California  Current  System. 

In  contrast,  coastal  modeling  studies  of  the  CCS  have  been  scarce.  McClain  et  al. 
(1986)  performed  the  first  limited  modeling  study  in  the  region,  using  ship-derived  winds 
to  produce  a  large  negative  wind  stress  curl  off  the  northwest  coast  of  the  IP,  resulting  in 
opposing  equatorward  and  poleward  surface  currents.  Recently  Batteen  et  al.  (1992) 
used  an  eddy  resolving  PE  model  with  both  uniform  and  variable  wind  stress,  the  latter 
computed  from  synoptic  surface  pressure  analyses,  to  produce  realistic  currents  as  well  as 
mesoscale  eddies.  This  study,  however,  covered  only  a  limited  region  of  the  CCS  (the 
northwest  coast  of  the  DP)  and  had  a  straight  coast  in  its  model  domain. 


The  objective  of  this  study  is  to  investigate  the  roles  of  seasonal  climatological 
wind  forcing  and  irregular  coastline  geometry  in  the  generation  of  currents,  eddies,  jets 
and  filaments  in  the  CCS.  The  high-resolution,  multi-level,  PE  model  of  Batteen  and 
Monroe  (1998)  will  be  used.  To  allow  larger  scale  eddies  and  elongated  filaments  to  be 
generated  and  the  model  to  reach  quasi-equilibrium,  the  model  will  be  allowed  to  run  for 
~3  years. 

This  study  is  organized  as  follows:  the  PE  model,  type  of  wind  forcing,  initial  and 
specific  experimental  conditions,  and  the  energy  analysis  technique  are  described  in 
chapter  n.  Results  of  the  model  simulations  and  comparisons  with  observations  are 
discussed  in  chapter  III.  A  summary  and  discussion  are  presented  in  chapter  IV. 


H.         MODEL  DESCRIPTION 


A.  MODEL  EQUATIONS 

The  PE  numerical  model  used  in  this  study  was  originally  a  coarse  resolution 
model  used  by  Haney  (1974)  to  study  closed  basins.  It  has  since  been  adapted  for  eddy- 
resolving,  limited  EBC  regions  with  open  northern,  western,  and  southern  boundaries  by 
Batteen  (1997)  The  multi-level,  non-adiabatic,  model  uses  beta-plane,  hydrostatic, 
Boussinesq,  and  rigid  lid  approximations  and  has  baroclinic  and  barotropic  velocity 
components.  Equations  governing  the  model  are  as  follows: 


-r  = -+fv-  AMV*u  +  KM—r 

dt      p0  dx     J         M              M  dz2 

(1) 

dv      - 1  dp                     ,              d2v 

dt      p0  dy     J         M              M  dz2 

(2) 

du     &     dw 
dx.     dy     dz 

(3) 

dp 

(4) 

p  =  p0[\-a(T-T0)  +  /?(S-SQ)] 

(5) 

dl ' =-ahv<t+k/J2 

dt          H              H  dz2 

(6) 

dt           "               H  dz2 

(7) 

In  the  above  equations,  /  is  time,  T  is  temperature,  S  is  salinity,  p  is  density,  and  p  is 
pressure.  A  right-handed  Cartesian  coordinate  system  (x,y,z)  is  used  where  x  points 
toward  shore,  y  alongshore,  and  z  upward.  The  corresponding  velocity  components  are 
(u,v,w).  Table  1  provides  a  list  of  symbols  found  in  the  model  equations,  as  well  as  values 
of  constants  used  throughout  the  study. 


A  space-staggered  B-scheme  (Arakawa  and  Lamb,  1977)  is  used  for  the  horizontal 
finite  differencing.  This  scheme  has  been  shown  by  Batteen  and  Han  (1981)  to  be 
appropriate  when  the  grid  spacing  is  approximately  on  the  same  order  as,  or  less  than,  the 
Rossby  radius  of  deformation.  The  horizontal  grid  spacing  is  9  km  in  the  east-west 
direction  and  11  km  in  the  north-south  direction,  while  the  internal  Rossby  radius  of 
deformation  is  -30  km. 

Consistent  with  Haney  (1974),  the  vertical  spacing  scheme  of  the  model 
concentrates  more  layers  above  the  thermocline  in  the  dynamically  active  portion  of  the 
ocean.  Specifically,  the  model  uses  ten  vertical  layers,  with  constant  z-levels,  at  depths  of 
10,  30,  75,  150,  250,  400,  600,  1226,  2283,  and  3656  m. 

The  model  domain  (Figure  1)  encompasses  the  west  coasts  of  the  IP  and  northwest 
Africa,  from  30°  N  to  42.5°  N  (1408  km  alongshore),  and  from  5°  W  to  17.5°  W  (1152 
km  cross-shore).  The  coastal  boundaries  of  the  model  domain  are  closed,  and  have  both 
the  tangential  and  normal  components  of  velocity  set  to  zero.  To  highlight  the  role  played 
by  climatological  wind  forcing,  bottom  topography  has  been  omitted  and  the  eastern 
boundary  is  modeled  as  a  vertical  wall.  The  constant  depth  used  in  the  model  is  4500  m. 

A  modified  version  of  the  radiation  boundary  conditions  of  Camerlengo  and 
O'Brien  (1980)  is  used  for  the  open  ocean  domain  boundaries  to  the  north,  south,  and 
west.  In  particular,  whether  a  boundary  grid  point  is  treated  as  an  inflow  point  or  an 
outflow  point  for  a  particular  prognostic  variable  is  determined  by  the  sign  of  a 
dynamically  computed  effective  group  velocity.  This  group  velocity  is  defined  as  the  ratio 
of  the  local  time  derivative  and  the  local  space  derivative  normal  to  the  boundary.  If  the 
boundary  grid  point  is  thereby  determined  to  be  an  inflow  point,  then  the  value  of  the 
prognostic  variable  is  set  to  its  value  at  the  previous  timestep.  If  the  boundary  point  is 
determined  to  be  an  outflow  point,  its  value  is  set  to  that  of  the  nearest  interior  point. 
Spatial  smoothing  is  also  applied  with  a  1-2-1  weighting  window  within  five  grid  points 
(-50  km)  of  the  open  boundaries. 


Biharmonic  lateral  heat  and  momentum  diffusion  is  used  in  the  model  with  the 
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same  choice  of  coefficients  (i.e.,  2.0  xlO      cm  /s)  as  in  Batteen  (1997).   Holland  (1978) 

showed  that  highly  selective  biharmonic  diffusion  acts  predominantly  on  submesoscales, 

while  Holland  and  Batteen  (1986)  found  that  baroclinic  mesoscale  processes  can  be 

damped  by  Laplacian  lateral  heat  diffusion.  As  a  result,  the  use  of  biharmonic  lateral 

diffusion  should  allow  mesoscale  eddy  generation  via  barotropic  (horizontal  shear)  and/or 

2 
baroclinic  (vertical  shear)  instability  mechanisms.  As  in  Batteen  (1997),  weak  (0.5  cm  /s) 

vertical  eddy  viscosities  and  conductivities  are  used.  Bottom  stress  is  parameterized  by  a 

simplified  quadratic  drag  law  (Weatherly,  1972),  as  in  Batteen  (1997). 

The  method  of  solution  is  straightforward  with  the  rigid  lid  and  flat  bottom 

assumptions    because    the    vertically    integrated    horizontal    velocity    is    subsequently 

nondivergent.   The  vertical  mean  flow  can  be  described  by  a  streamfunction  which  can  be 

predicted  from  the  vorticity  equation,  while  the  vertical  shear  currents  can  be  predicted 

after  the  vertical  mean  flow  is  subtracted  from  the  original  equations.   The  other  variables, 

i.e.  temperature,  salinity,  vertical  velocity,  and  pressure,  can  be  explicitly  obtained  from 

the  thermodynamic  energy  equation  (6),  salinity  equation  (7),  continuity  equation  (3),  and 

hydrostatic  equation  (4),  respectively  (For  more  complete  details  on  the  method  of 

solution,  see  the  Appendix). 

B.         TYPE  OF  WIND  FORCING 

Previous  experiments  by  Batteen  et  al.  (1992)  investigated  the  role  of  steady 
equatorward  winds  with  anticyclonic  wind  stress  curl  in  the  generation  of  features  off  the 
IP.  In  a  more  recent  study  Batteen  and  Monroe  (1998)  studied  the  contributions  of 
seasonal  wind  forcing  and  irregular  coastline  geometry  in  the  generation  of  eddies  and 
filaments  in  the  California  Current  system.  Following  Batteen  and  Monroe  (1998),  in  this 


study,  seasonal  wind  forcing  and  irregular  coastline  geometry  will  be  used  to  investigate 
the  generation  of  similar  features  in  the  CCS. 

To  explain  the  effects  of  seasonal  wind,  the  model  is  forced  from  rest  with 
climotalogical  wind  fields  from  a  2.5°  by  2.5°  grid  of  the  European  Centre  for  Medium 
Range  Weather  Forecasts  (ECMWF)  near-surface  wind  analyses  (Trenberth  et  al.,  1990). 
The  monthly  mean  stresses  based  on  twice  daily  wind  analyses  from  1980-1989  have  been 
interpolated  spatially  to  the  9  by  1 1  km  model  resolution  and  temporally  to  daily  wind 
values. 

Sample  wind  fields  used  can  be  seen  in  Figure  2,  which  shows  the  annual  migration 
of  the  Azores  Subtropical  High  from  the  south  in  the  winter  (e.g.,  Figure  2d),  to  its  most 
northern  extent  in  the  summer  (e.g.,  Figure  2b).  The  atmospheric  pressure  pattern  for 
November  (Figure  2d)  depicts  a  low  to  the  north  and  the  Azores  High  to  the  south,  which 
results  in  a  wind  divergence  at  -40°  N.  This  pattern  of  weakly  poleward  winds  north  of 
40°  N  and  equatorward  winds  to  the  south  continues  through  December.  During  January 
(not  shown)  and  February  (Figure  2a)  the  divergence  in  the  wind  field  migrates  poleward. 
By  March  (not  shown)  an  equatorward  component  in  the  wind  field  is  observed  along  the 
entire  domain.  The  strongest  equatorward  winds  are  discernible  from  July  (Figure  2b) 
through  August  (not  shown).  By  September  (Figure  2c)  the  winds  start  to  weaken 
throughout  the  domain,  and  divergence  in  the  wind  field  is  observed  in  the  north  in 
October  (not  shown).  In  November  (Figure  2d)  the  wind  divergence  returns  to  -40°  N. 

C.         EXPERIMENTAL  DESIGN 

The  design  of  the  model  experiment  is  as  follows.  The  model  is  forced  from  rest 
with  seasonal  ECMWF  winds.  The  initial  mean  stratification  used  are  annual 
climatological  temperature  fields  based  on  Levitus  and  Boyer  (1994)  centered  at  -37.5°  N 
(corresponding  to  the  center  of  the  model  domain).    The  temperatures  (°C)  used  for  the 


ten  levels  from  the  surface  to  4500  m  are  17.5,  17.3,  16.5,  15.0,  13.7,  12.5,  11.0,  8.52, 
3.59,  and  2.09,  respectively,  while  a  salinity  constant  of  34.7  is  used  for  all  levels. 

Like  all  major  EBC  systems,  the  CCS  is  a  region  of  net  annual  gain.  This  heat  gain 
occurs  because  of  relatively  low  cloud  cover  (compared  with  farther  offshore),  reduced 
latent  heat  flux,  and  downward  sensible  heat  flux  due  to  the  presence  of  cold  upwelled 
water  during  summer.  To  focus  the  experiment  on  wind  forcing  as  a  possible  mechanism 
for  the  generation  of  thermal  variability  in  the  CCS,  the  surface  thermal  forcing  in  the 
model  was  highly  simplified.  The  solar  radiation  at  the  sea  surface  S0,  was  specified  to  be 
the  summer  mean  and  CCS  mean  value.  On  the  other  hand,  the  sum  of  the  net  longwave 
radiation,  latent,  and  sensible  heat  fluxes,  Qb,  was  computed  during  the  model's 
experiments  from  standard  bulk  formulas  (Haney  et  al.  1978)  using  the  summer  and  CCS- 
mean  value  of  the  alongshore  wind  (above),  cloud  cover,  relative  humidity,  air 
temperature,  and  model-predicted  sea  surface  temperature.  The  initial  sea  surface 
temperature  was  chosen  so  that  the  total  heat  flux  across  the  sea  surface,  So  -  Qb,  was 
zero  at  the  initial  time.  Therefore,  the  only  surface  heat  flux  forcing  in  the  experiments 
was  that  which  developed  in  Qb  as  a  result  of  (wind  forced)  fluctuations  in  the  sea  surface 
temperature.  As  discussed  in  Haney  (1985),  such  a  surface  thermal  forcing  damps  the  sea 
surface  temperature  fluctuations  to  the  atmosphere  on  a  time  scale  of  the  order  of  100 
days.  Consequently,  sea  surface  temperature  fluctuations  that  develop  due  to  wind  forcing 
should  be  observed  long  before  they  are  damped  by  the  computed  surface  heat  flux. 
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m.       RESULTS  FROM  THE  MODEL  SIMULATION 


A.         SPIN-UP  PHASE 

On  day  one  of  the  model  (1  January),  the  Azores  High  is  stationed  near  its 
southernmost  position.  In  response  to  predominantly  equatorward  winds,  a  coastal 
equatorward  surface  current  develops  in  the  southern  end  of  the  domain  and  by  day  60 
extends  along  the  entire  coast  (not  shown). 

In  the  spring  as  the  Azores  High  migrates  north,  the  wind  intensifies  and 
transitions  to  equatorward  flow  over  the  entire  domain.  As  a  result,  increased  stress  is 
exerted  on  the  ocean  surface  creating  Ekman  transport  offshore  and  upwelling  of  cooler 
water  along  the  coast.  Upwelling  predominantly  occurs  in  the  south  where  stronger  winds 
exist  (e.g.,  Figure  3a),  but  by  mid-summer  is  evident  along  the  entire  coast  (e.g.,  Figures 
3b  and  4).  Upwelling  is  enhanced  near  capes  (e.g.,  Figure  3b),  where  both  the  alongshore 
component  of  the  wind  stress  and  the  coastal  current  velocity  are  at  local  maxima 
(Batteen,  1997). 

Below  the  equatorward  surface  current,  a  poleward  undercurrent  develops  first  at 
the  equatorward  side  and  then  at  the  poleward  side  of  the  model  domain.  The  coastal  jet 
extends  from  -200  m  depth  near  the  coast  to  -600-800  m  depth  offshore,  and  has  a  peak 
core  velocity  of  -30-45  cm/s  (Figures  5a  and  5b).  Along  the  west  coast  of  the  IP  (e.g., 
Figure  5  a)  and  Morocco  (not  shown)  the  undercurrent  develops  separate  cores  at  -400  m 
depth  and  1200  m  depth.  The  -400  m  depth  core  of  the  undercurrent  has  a  peak  velocity 
of -5-10  cm/s  and  remains  within  -20-30  km  of  the  coast.  The  -1200  m  depth  core  of  the 
undercurrent  has  a  peak  velocity  of -5-10  cm/s  and  generally  remains  within  -40-50  km  of 
the  coast. 
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In  the  Gulf  of  Cadiz  there  is  no  separation  of  the  undercurrent  cores.  The  core  of 
the  undercurrent  is  found  between  -400  and  1200  m  depth  and  has  speeds  of  5  cm/s 
(Figure  5b). 

As  the  coastal  jet  and  the  subsurface  undercurrent  become  fully  established  (for 
example,  see  Figures  5a  and  5b,  which  shows  the  structure  of  the  currents  prior  to 
meander  formation),  strong  vertical  and  horizontal  shears  develop  in  the  upper  layers 
between  the  opposing  currents.  This  shear  causes  the  flow  to  become  both  barotropically 
and  baroclinically  unstable  and  leads  to  the  formation  of  meanders  and  filaments. 
Meanders  occur  first  off  Cape  Ghir,  Cabo  da  Roca  and  Cabo  de  Sao  Vincent.  These 
meanders  intensify  and  develop  into  predominantly  cyclonic  cold  core  eddies,  which  in 
time  coalesce  with  other  eddies  to  form  larger  cyclonic  eddies  on  the  order  of  -150-250 
km.  By  day  285  (Figure  3c)  meanders  and  eddies  are  visible  throughout  the  coastal 
region. 

In  the  late  fall  as  the  Azores  High  migrates  south,  first  the  winds  and  then  the 
currents  begin  to  weaken.  In  time,  upwelling,  as  expected,  also  weakens.  By  the  end  of 
the  first  year  of  model  simulation  (Figure  3d),  the  equatorward  surface  flow  extends  as  far 
as    -350  km  offshore  and  has  taken  the  form  of  a  meandering  jet  embedded  with 
predominantly  cyclonic  eddies. 

B.         QUASI-EQUILBRIUM  PHASE 

Continuous  formation  and  sustainment  of  meanders,  eddies  and  filaments  indicate 
that  the  model  has  reached  a  quasi-equjlibrium  state  by  year  three.  Time-averaging  model 
output  fields  for  each  month  as  well  as  the  peak  upwelling  season  (July  through 
September)  in  year  three  provides  an  opportunity  to  see  the  complex  structure  and 
seasonal  variability  of  the  CCS. 

The  winter  model  results  (e.g.,  Figure  6a)  depict  a  predominantly  equatorward 
flow,  which  is  embedded  with  eddies  throughout  the  domain.    Along  the  IP,  a  surface 
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coastal  poleward  flow  with  speeds  as  high  as  -22.5  cm/s  is  visible  (Figures  6a  and  7a). 
Beneath  the  surface  (Figure  7a),  cores  of  poleward  flow  with  speeds  of  -10-17.5  cm/s  at 
-600  and  -1200  m  depth  are  discernible.  An  additional  third  core  of  poleward  flow  at 
-250  m  depth  with  a  speed  of -10  cm/s  is  also  formed  along  the  IP  in  the  vicinity  of -39° 
N  (Figure  7a). 

By  early  spring  an  equatorward  coastal  jet  appears  off  the  coast  of  Morocco  south 
of  32°  N  and  up  welling  resumes  (not  shown).  By  June  (Figure  6b)  upwelling  appears  all 
along  the  IP  and  a  strong  equatorward  current  with  speeds  of  -25-40  cm/s  meanders 
along  the  IP.  Off  Cabo  da  Roca  the  equatorward  current  weakens  as  it  flows  into  the  Gulf 
of  Cadiz.  South  of  35°  N  the  equatorward  current  becomes  a  coastal  jet  with  speeds  of 
-40  cm/s.  Downstream  of  Cabo  da  Roca  (Figure  6b)  and  Cape  Beddouzza  (not  shown) 
pairs  of  cyclonic  (Figures  6b  and  7b)  and  anticyclonic  (Figures  6b  and  7c)  eddies  of  O(150 
km)  in  size  and  extending  to  greater  than  1400  m  depth  are  discernible. 

In  the  summer  (Figure  6c)  upwelling  is  visible  along  the  western  and  southern 
coast  of  the  IP  as  well  as  along  the  coast  of  Morocco.  As  the  equatorward  jet  flows  along 
the  IP,  it  impinges  upon  Cabo  da  Roca  and  is  deflected  offshore  advecting  cold  water  with 
it,  resulting  in  the  formation  of  filaments  (e.g.,  Figure  6c)  off  Cabo  da  Roca  and  Cabo  de 
Sao  Vincent.  Filaments  are  also  discernible  off  Cabo  de  Santa  Maria  and  the  southern 
coast  of  Portugal  (Figure  6d).  The  filaments  typically  extend  from  -75  to  150  m  depth 
(not  shown).  Beneath  the  surface  (Figure  7d),  cores  of  poleward  flow  with  speeds  of -5- 
7.5  cm/s  at  -600  and  -1200  m  depth  are  evident.  An  additional  weaker  core  of  poleward 
flow  at  -250  m  depth  with  a  speed  of -5  cm/s  is  also  formed  off  the  west  coast  of  the  IP 
in  the  vicinity  of  37°  N. 

The  results  (Figure  8)  of  averaging  the  temperature  and  current  fields  for  the 
period  when  the  most  intense  upwelling  (i.e.,  July  through  September)  occurs  off  both  the 
IP  and  Morocco  shows  that  a  meandering  equatorward  current  offshore,  an  equatorward 
coastal  jet,  subsurface  poleward  flow  with  cores  at  -600  m  depth  and  -1200  m  depth, 
coastal  upwelling,  filaments  and  eddies  are  regular  features  of  the  CCS  during  the 
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upwelling  season.  Figure  8a  shows  the  surface  equatorward  current  offshore  of  the  ~ 
17.5°  C  isotherm  until  south  of  Cabo  da  Roca  where  the  16.5°  C  and  17.5°  C  isotherms 
separate.  South  of  Cabo  de  Sao  Vincent  the  offshore  equatorward  current  once  again 
splits  into  two  distinct  flows.  The  westernmost  flow  continues  equatorward  while  the 
easternmost  flow  meanders  cyclonically  and  joins  the  coastal  jet  in  the  Gulf  of  Cadiz.  Off 
Cape  Beddouzza  the  offshore  flow  and  coastal  current  rejoin,  and  continue  down  the 
coast  of  Morocco.  The  offshore  current  has  speeds  of -15-45  cm/s  and  extends  to  -1000 
m  depth  (Figures  8b  and  8c).  The  coastal  jet  has  peak  speeds  of  -45  cm/s  and  extends  to 
-250  m  depth  near  the  coast  (Figure  8c).  Upwelling  is  visible  all  along  the  IP  coast  as 
well  as  along  the  coast  of  Morocco  (Figure  8a).  Cold,  upwelled  filaments  are  also 
discernible  downstream  of  prominent  capes  (Figure  8a).  The  filaments  typically  extend 
-125  km  offshore,  have  widths  of  -75  km  and  extend  to  -100  m  depth  (not  shown). 
Figure  8b  shows  the  often  fuzed  -250  m  depth  and  600  m  depth  undercurrent  cores  found 
between  Cabo  da  Roca  and  Cabo  de  Sao  Vincent  as  well  as  the  -1200  m  depth 
undercurrent  core.  The  -250  m  and  600  m  depth  cores  reach  speeds  of -10  cm/s  while 
the  1200m  depth  core  reaches  speeds  of -2.5-5  cm/s.  Figure  8d  shows  the  -1200  m  depth 
undercurrent  as  it  flows  along  the  entire  coast  from  Morocco  into  the  Gulf  of  Cadiz  and 
along  the  IP  during  the  peak  upwelling  season.  Offshore  of  the  undercurrent  are  both 
cyclonic  and  anticy clonic  eddies  of  O(150)  km  in  size  that  extend  to  at  least  1000  m  depth 
(Figures  8d  and  8e). 

Figures  9a  and  9b  show  horizontal  maps  of  upper  layer  mean  kinetic  energy 
(MKE)  and  eddy  kinetic  energy  (EKE)  averaged  during  the  peak  upwelling  season. 
Holland  et  al.  (1983)  showed  that  maps  of  MKE  and  EKE  can  be  used  to  locate  sources 
of  mean  and  eddy  energy.  From  Figure  9a  it  can  be  seen  that  high  values  of  MKE  follow 
the  equatorward  jet  and  that  the  highest  values  of  MKE  appear  near  areas  of  tightened 
temperature  gradients  off  Cabo  da  Roca,  in  the  Gulf  of  Cadiz  and  off  Cape  Beddouzza. 

A  comparison  of  Figures  8a  and  9b  shows  maximum  values  of  EKE  in  the  vicinity 
of  the  large  cyclonic  meander  that  extends  downstream  from  Cabo  da  Roca.   High  values 
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are  also  formed  along  the  axes  of  the  equatorward  jet  on  the  west  coast  of  the  Gulf  of 
Cadiz  and  offshore  of  Cape  Ghir  and  Cape  Beddouzza.  The  lowest  coastal  values  of  EKE 
are  along  the  southern  coast  of  the  IP,  indicating  that  this  region  is  not  a  primary  source  of 
eddy  energy. 

A  comparison  of  Figures  8a,  9a  and  9b  shows  that  overall,  the  regions  of 
maximum  MKE  and  EKE  are  found  along  the  axes  of  the  equatorward  jet.  This  is 
consistent  with  the  results  of  the  model  simulation,  which  showed  that  the  eddies  are 
generated  from  instabilities  of  the  mean  equatorward  current  and  subsurface  undercurrent 
via  baroclinic  and/or  barotropic  instability  processes. 

By  Fall  (e.g.,  Figure  6e)  upwelling  has  weakened  and  the  coastal  jet  is  no  longer 
visible  north  of  Cabo  da  Roca.  Offshore,  a  meandering  equatorward  current  embedded 
with  westward  propagating,  cyclonic  eddies  is  discernible.  Farther  south,  the  coastal  jet 
leaves  the  coast  at  Cape  Beddouzza  and  is  entrained  in  a  large  cyclonic  eddy  0(200  km) 
prior  to  continuing  south  along  the  coast  of  Morocco.  By  winter,  poleward  surface  flow 
with  speeds  of -15-25  cm/s  has  developed  north  of  41°  N.  The  flow  is  found  within  ~60 
km  of  the  coast  and  extends  to  -200  m  depth  (Figure  7e). 

C.         COMPARISONS  OF  MODEL  RESULTS  WITH  OBSERVATIONS 

To  ensure  model  results  are  consistent  with  realistic  features  of  the  CCS,  the 
model  results  are  compared  to  available  field  and  satellite  observations.  Since  this  is  an 
idealized  process-oriented  study  and  not  a  model  hindcast,  direct  comparisons  are  not 
valid;  however,  it  can  be  shown  that  the  model's  phenomenological  behavior  is 
qualitatively  similar  to  observed  data.  The  discussion  in  the  ensuing  paragraphs  first 
details  general  CCS  features  and  then  specific  magnitudes  and  characteristics  of  features 
generated  along  the  western  coasts  of  the  IP  and  Northwest  Africa.  Overall,  there  exists 
impressive  similarities  between  the  results  of  the  model  simulation  and  the  major 
characteristics  of  the  CCS. 


15 


1.  Comparison  of  Ocean  Currents 

The  major  characteristics  of  the  modeled  Canary  Current  off  western  Iberia  and 
northwestern  Africa  are  consistent  with  the  few  field  observations  and  studies  available. 
Satellite  imagery  also  reveals  a  system  of  meanders,  eddies,  and  along-shelf  jets,  on  a 
range  of  length  scales  that  correspond  closely  with  model  output. 

a.  Equatorward  Flow 

Using  climatological  wind  data  and  a  realistic  coastline,  the  model  develops 
an  ocean  current  with  features  similar  to  those  observed  in  the  CCS.  The  broad, 
equatorward  flow  found  throughout  the  model  domain  matches  the  general  description  of 
classical  EBCs  by  Wooster  and  Reid  (1963).  The  sustained  magnitude  of  the  model's 
equatorward  flow  also  compares  favorably  with  available  data.  For  example,  DMAH/TC 
(1988)  shows  a  mean  predominantly  equatorward  flow  of -0.5  to  0.7  kts  (-25-35  cm/s), 
which  is  consistent  with  velocities  found  in  the  model  simulation  (e.g.,  Figures  8a  and  8c). 

b.  Poleward  Flow 

The  model  establishes  a  poleward  undercurrent  with  cores  at  -600  m  depth 
and  -1200  m  depth  throughout  the  domain  (Figures  5a,  7a,  8b,  and  8c).  In  addition,  a 
third  core  of  poleward  flow  at  -250  m  depth  is  seen  near  37°  N  along  the  west  coast  of 
the  IP  predominantly  during  June,  July  and  August  (e.g.,  Figure  7d)  and  near  39°  N  during 
November  and  January  (e.g.,  Figure  7a).  The  depths  and  magnitudes  of  the  600  m  and 
1200  m  depth  model's  poleward  flow  cores  are  consistent  with  Ambar  (1980),  which 
shows  these  cores  located  between  600-900  m  and  1100-1200  m  depth  along  the  entire 
western  coast  of  the  IP.  The  depth  and  speed  of  the  250  m  model's  poleward  flow  is 
consistent  with  Haynes  and  Barton  (1990),  which  shows  this  core  located  between  100 
and  600  m  depth  off  the  IP. 
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Recent  studies  by  Frouin  et  al.  (1990)  and  Haynes  and  Barton  (1990), 
using  high  resolution  infra-red  images  from  NOAA7  and  NOAA9  and  numerous  in-situ 
measurements,  reveal  the  existence  of  surface  poleward  flow  off  the  northern  coast  of 
Portugal  and  northwestern  Spain.  These  studies  show  a  current  -25-40  km  wide,  that 
extends  to  -200  m  depth  and  reaches  speeds  of -20-25  cm/s.  Similarly,  in  response  to  the 
shift  in  the  climatological  winds  from  northerly  to  southwesterly  late  in  the  year  (see 
Figure  2),  the  model  generates  surface  poleward  flow  in  the  northern  end  of  the  domain. 
As  seen  in  Figures  6a,  7a  and  7e  surface  poleward  flow  with  speeds  of  15-25  cm/s 
develops  within  60  km  of  the  coast  and  to  a  depth  of -200  m. 

2.  Comparison  of  Eddies 

High  resolution  satellite  imagery  and  observations  by  Fiuza  (1984)  show  pairs  of 
mesoscale  eddies  on  the  order  of  100  km  off  the  west  coast  of  the  IP  (Figure  10).  This 
study  confirms  the  model's  generation  of  numerous  eddies  on  the  order  of  100  km  in  the 
CCS. 

Preferred  locations  of  instability  onset  (eddy  generation)  occur  where  the  coastline 
significantly  changes,  i.e.,  near  Cabo  da  Roca,  Cabo  de  Sao  Vincent,  the  Gulf  of  Cadiz, 
and  Cape  Beddouzza/Cape  Ghir.  A  comparison  of  the  results  of  Fiuza  (1984)  (Figure  10), 
with  the  results  of  the  model  simulation  (e.g.,  Figures  3c,  3d,  6  and  11)  shows  that  the 
model  generates  eddies  in  similar  locations. 

3.  Comparison  of  Upwelling 

Upwelling  varies  by  season  throughout  the  domain  as  a  result  of  the  seasonally 
oscillating  winds  seen  in  Figure  2.  Observations  show  that  upwelling  appears  first  off 
Cape  Ghir/Cape  Beddouzza  in  early  spring  as  seen  in  Van  Camp  et  al.  (1991),  and  then 
progressively  to  the  north  as  the  Azores  High  moves  up  the  Moroccan  coast.  This 
progression  of  upwelling  is  well  simulated  in  the  model  (e.g.,  Figures  11a  and  lib). 
Observed  upwelling  off  Portugal  from  -43°  N  to  Cabo  de  Sao  Vincent  in  the  south 
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(Figure  12a)  and  along  the  southern  coast  of  Portugal  to  Cabo  de  Santa  Maria,  -8°  W 
(Folkard  et  al.,  1997)  is  also  well  depicted  by  the  model  (e.g.,  Figure  lie).  Observations 
of  peak  upwelling  during  July,  August  and  September  are  shown  in  Figures  12b  and  12c 
while  Figures  8a,  lid  and  lie  show  the  model  depicting  the  most  intense  upwelling  during 
this  period  as  well. 

4.  Comparison  of  Filaments 

Several  studies  by  Fiuza  (1983),  Barton  (1986),  Fiuza  and  Sousa  (1989),  Folkard 
et  al,  (1997),  Van  Camp  et  al.  (1991)  and  Hagen  et  al.  (1996)  have  confirmed  the 
existence  of  cold  filaments  extending  off  the  west  coasts  of  the  EP  and  Morocco  depicted 
by  the  model. 

a.  West  Coast  of  Iberia 

Fiuza  (1983)  was  the  first  study  to  indicate  the  existence  of  filaments  off 
the  west  coast  of  the  IP.  Subsequent  studies  by  Fiuza  and  Sousa  (1989)  using  CTD  data 
and  satellite  imagery  confirmed  the  existence  of  filaments  including  a  large  filament  off 
Cabo  de  Roca.  In  addition,  recent  SST  observations  show  several  good  images  of 
filaments  occurring  along  the  entire  west  coast  of  the  IP  (e.g.,  Figure  12b)  with  a 
particularly  long  filament  at  40.5°  N.  The  model  depicts  filaments  in  many  of  the  same 
locations  as  well  (Figures  lid,  lie,  1  If  and  1  lg).  A  filament  exists  near  -41°  N  (Figures 
lid  and  1  le)  while  a  rather  strong  filament  is  also  evident  at  -39°  N  (Figures  1  Id,  1  If  and 
1  lg)  and  downstream  of  Cabo  da  Roca  (Figures  1  Id  and  1  le). 

b.  South  Coast  of  Iberia/Gulf  of  Cadiz 

While  investigating  the  waters  of  southern  Iberia  in  the  Gulf  of  Cadiz, 
Folkard  et  al.  (1997)  observed  cold,  upwelled  filaments  extending  southward  from  the 
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coast.  The  model  results  shown  in  Figures  6d,  6e,  1  If  and  llg  depict  these  filaments  as 
well. 

c.  Cape  Ghir/Northwest  Africa 

Satellite  observations  from  Van  Camp  et  al.  (1991)  and  Hagen  et  al. 
(1996)  show  a  filament  of  cold,  upwelled  water  extending  off  Cape  Ghir  near  31°  N  for 
about  200  km.  Consistent  with  these  observations,  the  model  depicts  this  filament  as  an 
equatorward  meandering  coastal  jet  which  transports  cold,  upwelled  water  offshore 
(Figure  lie).  Observations  show  that  the  coastal  jet  and  subsequent  frontal  zone  caused 
by  the  advected  cold  water  results  in  the  generation  near  -12°  W  of  a  cyclonic  eddy-like 
feature  on  the  southern  side  of  the  jet,  while  an  anticyclonic  eddy  forms  on  the  northern 
side.  The  model  results  (Figure  lie)  show  that,  consistent  with  the  observations,  on  the 
northern  side  of  the  jet  is  an  anticyclonic  eddy  centered  at  -12°  W,  33°  N  and  on  the 
southern  side  of  the  jet  a  cyclonic  eddy  at  -12°  W,  31.5°  N. 

In  summary,  the  variety  of  length  scales  and  life  time  scales,  in  conjunction 
with  the  differences  in  field  observations  themselves,  emphasize  the  complex  flow  region 
of  the  CCS.  Despite  this  complexity,  the  numerical  model  results  show  striking  similarities 
to  available  field  observations  of  surface  and  subsurface  currents,  upwelling,  filaments, 
meanders,  and  cyclonic  and  anticyclonic  eddies. 
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IV.       SUMMARY 

The  objective  of  this  study  was  to  investigate  the  role  of  seasonal  climatological 
wind  forcing  in  the  generation  of  currents,  eddies,  jets,  and  filaments  in  the  CCS.  Toward 
this  end,  a  high  resolution,  multi-level,  PE  model  using  a  realistic  coastline  was  forced 
from  rest  using  spatially  and  temporally  varying  winds.  The  cyclic  migration  of  the  Azores 
High,  and  the  subsequent  role  it  has  on  the  seasonal  variability  of  alongshore  winds,  was 
shown  to  play  a  significant  role  in  the  generation,  evolution  and  sustainment  of  mesoscale 
features  in  the  CCS. 

An  equatorward  surface  current  and  a  subsurface  poleward  current  with  cores  at 
-600  m  and  1200  m  depth  developed  and  were  maintained  throughout  the  year.  A  third 
core  of  poleward  flow  was  also  generated  at  -250  m  depth  at  -37°  N  from  June  through 
August  and  at  -39°  N  during  January  and  November.  Along  the  northern  IP,  surface 
poleward  flow  developed  in  the  late  fall  and  lasted  until  late  spring.  The  opposing  nature 
of  the  equatorward  and  subsurface  poleward  currents  resulted  in  both  horizontal 
(barotropic  instability)  and  vertical  (baroclinic  instability)  shear.  This  shear  caused  the 
currents  to  become  unstable,  meander,  and  eventually  lead  to  the  formation  of  cyclonic 
and  anticyclonic  eddies.  As  meanders  intensified  along  the  coast,  cold,  upwelled  water 
was  advected  offshore  forming  filaments. 

Model  results  also  depicted  a  seasonal  upwelling  cycle  tied  to  the  migrating  Azores 
High.  Upwelling  began  in  the  spring  along  the  southern  coast  of  Morocco,  progressed 
northward,  and  by  mid-summer  was  visible  along  the  coast  throughout  the  domain. 
During  the  upwelling  season  a  surface  equatorward  coastal  jet  overlying  the  poleward 
undercurrent  was  generated.  As  the  coastal  jet  reached  its  maximum  velocity  during  the 
peak  upwelling  season  (July  through  September),  eddy  formation  along  the  coast 
increased.  Seasonal  filament  generation  also  occurred  during  and  just  after  the  upwelling 
season. 
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Coastline  irregularities  also  appeared  to  play  a  significant  role  in  the  formation  of 
mesoscale  features.  Specifically,  the  capes  acted  as  barriers  that  forced  the  coastal 
equatorward  flow  offshore.  This  continuous  offshore  deflection  caused  some  meanders 
and  eddies  to  appear  as  quasi-permanent  rather  than  seasonal  features.  Capes  were  also 
areas  of  enhanced  upwelling  and  were  significant  in  the  location  and  formation  of 
filaments. 

A  comparison  of  the  currents,  eddies,  jets,  and  filaments  generated  by  the  model 
with  available  observations  showed  that  the  model  successfully  simulated  the  location, 
size,  and  velocity  of  these  features.  Of  the  features  generated,  the  poleward  undercurrent 
and  its  cores  are  of  particular  interest.  Until  now,  it  has  generally  been  accepted,  that 
these  cores  were  the  result  of  Mediterranean  Outflow  influence  (e.g.,  Ambar,  1980).  This 
study  showed,  however,  that  the  cores  could  be  generated  by  wind  forcing  alone. 

Future  modeling  efforts  to  enhance  the  simulation  and  resolution  of  CCS  features 
should  include  the  addition  of  thermohaline  gradients,  bottom  topography  and  MO. 
Observations  have  shown  that,  depending  on  prevailing  wind  direction,  an  inflow  or 
outflow  through  the  Strait  of  Gibraltar  can  exist,  which  can  affect  the  flow  in  the  Gulf  of 
Cadiz  (e.g.,  Folkard  et  al.,  1997).  Opening  the  Strait  of  Gibraltar  into  the  Mediterranean 
Sea  would  also  allow  the  study  of  density-driven  undercurrent  components.  With  the 
addition  of  bottom  topography,  the  topographical  steering  of  MO  and  anchoring  of  coastal 
features  could  be  studied. 

Using  ECMWF  winds  of  relatively  coarse  spatial  resolution,  this  study  focused  on 
the  role  of  seasonal  winds  in  the  CCS.  It  would  be  of  interest  to  investigate  this  role  using 
winds  of  higher  spatial  and  temporal  resolution.  The  use  of  winds  with  spatial  resolution 
down  to  ~1  km  and  temporal  resolution  of  ~1  day  would  allow  the  oceanic  response  to 
wind  events,  relaxations,  and  reversals  to  be  studied. 

The  lack  of  readily  available  data  along  the  Moroccan  and  IP  coasts  was  a  concern 
to  this  study.  Although  several  studies  discuss  briefly  the  surface  temperature,  salinity, 
and  density  patterns  off  the  EP,  and  northwest  Africa  and  for  the  Cape  Ghir  Filament,  little 
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discussion  is  available  about  the  basics.  Few  observations  concerning  the  velocity  and 
structure  of  equatorward  flow  exist.  Fewer  still  are  observations  of  surface  poleward  flow 
along  the  IP  and  subsurface  observations  of  filaments,  meanders  and  eddies.  In  short,  this 
area  remains  a  prime  candidate  for  a  refreshing  update  of  the  generalization  given  the 
Canary  Current  by  Wooster  and  Reid  (1963). 

Nevertheless,  the  results  from  this  experiment  show  surprising  similarities  to  the 
large-scale  observational  data  available  in  the  area  studied.  The  results  also  support  the 
hypothesis  that  wind  forcing  and  coastal  irregularities  play  a  key  role  in  the  generation, 
evolution,  and  sustainment  of  the  currents,  meanders,  eddies,  and  filaments  of  the  CCS. 
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Figure  1.  The  model  domain  for  the  Canary  Current  System  (CCS)  is  bounded  by  30°  N 
to  42.5°  N,  5°  W  to  17.5°  W.  The  cross-shore  (alongshore)  resolution  is  9  km  (11  km). 
Geographic  locations  and  prominent  features  are  labeled. 
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Figure  2.   Climatological  (1980-1989)  ECMWF  winds  in  m/s  for:  (a)  February,  (b)  July, 
(c)  September,  and  (d)  November.  Maximum  wind  vector  is  10  m/s. 
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Figure  3.  Temperature  contours  and  velocity  vectors  at  30  m  depth  at  days  (a)  150,  (b) 
195,  (c)  285  and  (d)  360.  The  contour  interval  is  1°C.  To  avoid  clutter,  the  velocity 
vectors  are  plotted  every  third  grid  point  in  the  cross-shore  direction  and  every  fourth  grid 
point  in  the  along-shore  direction.  Maximum  current  velocity  is  100  cm/s. 
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Figure  4.  Cross-shore  sections  of  isotherms  on  day  195  at  (a)  33°  N  and  (b)  39°  N. 
Contour  interval  is  1°  C. 
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Figure  5.  Cross-shore  sections  of  meridional  velocity  (v)  on  day  240  at  (a)  39.5  °  N,  and 
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Figure  6.  Temperature  contours  (b-e)  and  velocity  vectors  (a-e)  at  30  m  depth  in  the  third 
year  of  the  model  simulation  time-averaged  over  the  months  of  (a)  January,  (b)  June,  (c) 
August,  (d)  September  and  (e)  November.  The  contour  interval  is  1°  C.  Maximum 
current  velocity  is  100  cm/s. 
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Figure  7.  Cross-shore  sections  of  meridional  velocity  (v)  in  the  third  year  of  model 
simulation  time-averaged  for  (a)  January  at  39.4°  N,  (b)  June  at  38°  N  and  (c)  38.6°  N,  (d) 
August  at  37.5°  N  and  (e)  December  at  41.5  °  N.  Solid  lines  indicate  poleward  flow,  while 
dashed  lines  indicate  equatorward  flow.  The  contour  interval  is  2.5  cm/s  for  poleward 
flow  in  (a),  (d)  and  (e),  5.0  cm/s  for  poleward  flow  in  (b)  and  (c),  and  5  cm/s  for 
equatorward  flow. 
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Figure  8.  Time-averaged  plots  for  the  peak  upwelling  season  (July-September)  of 
temperature  contours  and  velocity  vectors  (a)  at  30  m  depth  and  (d)  at  1226  m  depth,  and 
a  cross-shore  section  of  velocity  (v)  at  (b)  37.5°  N,  (c)  39°  N  and  (e)  37.4°  N  .  The 
contour  interval  is  1°  C  in  (a).  The  maximum  current  velocity  is  100  cm/s  in  (a)  and  40 
cm/s  in  (d).  In  (b),  (c)  and  (e)  solid  lines  indicate  poleward  flow,  while  dashed  lines 
indicate  equatorward  flow.  The  contour  interval  for  poleward  flow  is  2.5  cm/s  in  (b)  and 
(c)  and  5  cm/s  in  (e).  The  contour  interval  for  equatorward  flow  is  5  cm/s. 
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Figure  9.  (a)  Mean  and  (b)  eddy  kinetic  energy  at  30  m  depth  for  the  peak  upwelling 
season  (July  through  September).  Contour  interval  is  100  cm2/s2. 
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Figure  10.     Computer  representation  of  AVHRR  NOAA7  Sea  Surface  Temperatures 
(SST)  with  observed  eddy  formations  from  5  August  1982.  "After  Fiuza  1984" 
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Figure  11  Temperature  contours  and  velocity  vectors  at  30  m  showing  upwelling  in  the 
third  year  at  day  (a)  137,  (b)  152,  (c)  173,  (d)  209,  (e)  263,  (f)  290,  and  (g)  317.  The 
contour  interval  is  1°  C.  The  maximum  current  velocity  is  100  cm/s. 
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Figure  12.  AVHRR  sea  surface  temperatures  (SST)  from  (a)  21  June  1995  showing 
increased  cold  SST  southwest  of  Cabo  de  Sao  Vincent  caused  by  intensified  up  welling  and 
advection  of  upwelled  water  around  the  cape  and  (b)  27  July  1993  showing  strong 
upwelling  and  numerous  filaments  off  the  west  coast  of  the  IP.  In  (c)  SST  show  intense 
upwelling  north  of  Cabo  da  Roca  on  29  July  1977  Figures  (a)  and  (b)  are  from  the 
Remote  Sensing  Data  Analysis  Service  (RSDAS),  Plymouth  Marine  Laboratory,  UK. 
Figure  (c)  is  from  Fiuza  (1982). 
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APPENDIX.  METHOD  OF  SOLUTION 

Equations  (1)  through  (7)  comprise  a  closed  system  of  seven  scalar  equations  and 
seven  unknowns,  u,  v,  w,  p,  p,  T,  and  S.  The  variables,  w,  v,  T,  and  S  are  prognostic 
variables  whose  time  rates  of  change  are  predicted  from  (1),  (3),  (6)  and  (7),  respectively. 
Although  the  diagnostic  variables  w,  p,  and  p  can  be  determined  from  (3),  (4),  and  (5), 
respectively,  additional  constraints  are  imposed  on  p  and  w  by  the  choice  of  the  rigid  lid 
boundary  conditions.  The  vertically  integrated  pressure  can  no  longer  be  obtained  by 
integrating  the  hydrostatic  equation  (4)  for  the  free  surface,  and  the  vertically-integrated 
horizontal  velocity  is  subsequently  constrained  to  be  non-divergent,  i.e., 


du     dv 


Li***)*-0-  (A1) 


which  is  obtained  by  integrating  (3)  and  applying  the  vertical  boundary  conditions  where 
s  is  a  dummy  variable  representing  the  vertical  coordinate. 

For  any  quantity  q,  let  its  vertical  average  be  denoted  by  q  and  its  departure 
(vertical  shear)  by  q' .  From  (Al)  the  vertical  mean  flow  can  then  be  described  by  a 
streamfunction  if/ ,  such  that: 

1  dw 
»  =  -H^    •  (A2> 

v  =  -L^         .  (A3) 

H  dx.  ' 

The  streamfunction  y/  is  predicted  from  the  vorticity  equation,  which  is  derived  by 

applying  the  curl  operator  to  the  vertical  average  of  (1)  and  (2),  and  then  using  (A2)  and 
(A3),  the  vorticity  equation  becomes 
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(A4) 
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dx\Hl-H     J     dy\H}-H    ) 


where  G  and  F  represent  the  collected  contributions  of  the  nonlinear  and  viscous  terms 
from  equations  (1)  and  (2). 

The  vorticity  equation  (A4)  is  solved  by  obtaining  an  updated  value  of  £  by 

application  of  the  leapfrog  (or  every  11  time  steps,  the  Euler-backward)  time-differencing 
scheme.  The  associated  value  of  y/  can  then  be  obtained  from: 


f-if^l+: 


HVdx- 


dy/  dHl      dy/  dH 

+  — r —  + 


H  \  dy    )      dx     dx         dy     dy 


(A5) 


which  is  an  elliptic  equation.  A  solution  to  (A5)  is  fully  prescribed  by  specifying  the  values 
of  if/  on  the  open  and  closed  boundaries  of  the  model  domain.  Currently,  to  solve  (A5), 

the  model  uses  an  elliptic  solver  when  there  are  no  variations  in  coastline  geometry  and/or 
topography,  and  successive  over-relaxation  techniques  when  there  are  variations  in 
coastline  geometry  and/or  topography. 

The  vertical  shear  current  (it',  v')  is  predicted  from  (1)  and  (2)  after  subtracting 
the  vertical  mean  flow.  The  results  are: 

(A6) 


a    pQ  dx 


dv'      -I  dp' 
a       Pq   dy 


d2v' 


PoH 


-     r 
+  G-G 


dz  p0H  (A7) 

In  (A6)  and  (A7),  p',  which  represents  the  departure  of  the  pressure  from  the  vertical 
average,  is,  using  (4),  expressed  in  terms  of  p  as: 


o  fo 


u  i     v    r  v 

p'=jpgds-  —  j    \pgds 

z  "  -H\z 


(A8) 


dz. 


H 
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The  method  of  solution  consists  of  predicting  V2y/,y/,u' ,v' ,T,  and  S  from  (A4), 

(A5),  (A6),  (A7),    (6)  and  (7),  respectively.  The  total  current  is  then  obtained  by  adding 
the  vertical  shear  part  to  the  vertical  average  part,  after  the  latter  is  obtained  from  y/ 

using  (A2)  and  (A3).  The  diagnostics  p,  w,  and  p'  are  then  obtained  explicitly  from  the 
equation  of  state  (5),  continuity  equation  (Al),  and  hydrostatic  relation  (A8),  respectively. 
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Table  1.  Values  of  Constants  Used  in  the  Model 


Constant 

Value 

Definition 

To 

278.2°K 

Constant  Reference  Temperature 

So 

34.7 

Constant  Reference  Salinity 

A> 

1.0276  gm  cm3 

Density  of  Sea  Water  At  T0  and  S0 

a 

2.4  x  lCr4  (°K)"J 

Thermal  Expansion  Coefficient 

fi 

7.5  x  1CT4 

Saline  Expansion  Coefficient 

K 

10 

Number  of  Levels  In  Vertical 

Ax 

9.0  x  105cm 

Cross- Shore  Grid  Spacing 

Ay 

1.1  x  106cm 

Alongshore  Grid  Spacing 

H 

4.5  x  105cm 

Total  Ocean  Depth 

At 

800  s 

Time  Step 

/o 

0.86  x  104  s"1 

Mean  Coriolis  Parameter 

8 

980  cm  s2 

Acceleration  of  Gravity 

AM 

2  x  1017  cm4  s"1 

Biharmonic  Momentum  Diffusion  Coefficient 

Ah 

2  x  1017  cm4  s"1 

Biharmonic  Heat  Diffusion  Coefficient 

KM 

0.5  cm2  s"1 

Vertical  Eddy  Viscosity 

KH 

0.5  cm2  s"1 

Vertical  Eddy  Conductivity 
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